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Abstract 

The complete gauge-invariant set of the one-loop QED corrections to the parity-nonconserving (PNC) 
amplitude in cesium and francium is evaluated to all orders in aZ using a local form of the Dirac-Fock 
potential. The calculations are performed in both length and velocity gauges for the absorbed photon and 
the total binding QED correction is found to be —0.27(3)% for Cs and —0.28(5)% for Fr. Moreover, a high- 
precision calculation of the electron-correlation and Breit-interaction effects on the 7s-8s PNC amplitude 
in francium using a large-scale configuration-interaction Dirac-Fock method is performed. The obtained 
results are employed to improve the theoretical predictions for the PNC transition amplitude in Cs and Fr. 
Using an average value from two most accurate measurements of the vector transition polarizability, the 
weak charge of ^^'^Cs is derived to amount to Qw = — 72.65(29)exp(36)thcor- This value deviates by 1.1a 
from the prediction of the standard model. The values of the 7s-8s PNC amplitude in ^^Spj. ^^^^ siOpj. 
obtained to be -15.49(15) and -14.16(14), respectively, in units of ixlO^^^ {-Qw)/N a.u. 

PACS numbers: 11.30.Er, 31.30.Jv, 32.80. Ys 
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I. INTRODUCTION 



Measurements of the parity nonconservation (PNC) effects in atoms provide sensitive tests of 
the standard model (SM) and impose constraints on physics beyond it The 6s-7s PNC 

amplitude in ^^'^Cs jsl] remains one of the most effective tool for such investigations. The mea- 
surement of this amplitude to a 0.3% accuracy lyilsl has stimulated a reanalysis of the theoretical 



flSH]. First, it was found flHIli 



that the role of the Breit 



predictions given in Refs. , - 

TT 

interaction had been underestimated. Then, it was pointed out II14I1 that the QED corrections may 
be comparable with the Breit corrections. The numerical evaluation of the vacuum-polarization 
(VP) correction lisl led to a 0.4% increase of the 6s-7s PNC amplitude in ^^^Cs, which resulted 
in a 2.2a deviation of the weak charge of ^•^^Cs from the SM prediction. This has triggered a great 
interest to calculations of the complete one-loop QED corrections to the PNC amplitude. 

While the VP contribution can easily be evaluated to a high accuracy within the Uehling ap- 
proximation, the calculation of the self-energy (SB) contribution is a much more demanding prob- 
lem (here and below we imply that the SE term embraces all one-loop vertex diagrams as well). 



In the plane wave approximation, that corresponds to zeroth order in aZ, it was derived in Refs. 

This correction, whose relative value equals to —a/ {2n), is commonly included in the 



In the pla 



definition of the nuclear weak charge. The «Z-dependent part of the SE correction to the PNC 



matrix element between s and p states was evaluated in Refs. yj 



1911 . These calculations, which 



are exact to first order in aZ and partiallv include higher-order binding effects, yield the binding 
SE correction of -0.9(1)% 013] and -0.85% The corresponding total binding QED cor- 
rection was found to amount to —0.5% and —0.43%, respectively. Despite this restored agreement 
with SM, the status of the QED correction to PNC in ^'^^Cs could not be considered as resolved 
until a complete aZ-dependence calculation of the SE correction to the 6s-ls transition amplitude 
is accomplished. This is due to the following reasons. First, in case of cesium {Z = 55) the pa- 
rameter aZ ^ 0.4 is not small and, therefore, the higher-order corrections, which are bey ond the 
A{aZ)lia{Xc/ Rnnc) term liigj . can be significant. Second, because the calculations fisl fiol l2(ll 
are performed for the PNC matrix element only, they do not include other SE diagrams which 
contribute to the 6s-7s transition amplitude. For instance, these calculations do not account for di- 
agrams in which the virtual photon embraces both the weak interaction and the absorbed photon. 
Third, strictly speaking, the PNC matrix element between the states of different energies is not 
gauge invariant. Despite the gauge-dependent part is suppressed by the small energy difference 
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Hi, 



11911 . estimates of the uncertainty in the definition of the PNC diagrams may fail due to unphysical 
origin of the gauge-dependent terms. 

The firs^step towards a complete aZ-dependence calculation of the SE correction was done 
in Ref. [21], where the SE correction to the 2s-2pi/2 PNC matrix element in H-like ions was 
evaluated. This matrix element was chosen to deal with the simplified gauge-invariant amplitude. 
The results of that work agree with those of Refs. However, as was stressed there, no 

claims can be made about the applicability of these results to the 6s-ls PNC transition in neutral 
cesium. 

Finally, the whole gauge-invariant set of the one-loop QED corrections to the 6s-ls PNC transi- 
tion amplitude in cesium was evaluated in Ref. 1221 . This calculation showed that the contributions 
of all SE diagrams are of the same order of magnitude (in both length and velocity gauges) and 
the final result arises through a delicate cancellation of individual terms, none of which can be 
neglected. The binding SE correction was obtained to amount to —0.67(3)% whereas the total 
binding QED correction is —0.27(3)%. 

Recently, the one-loop radiative corrections to the 6s-ls PNC amplitude in cesium were reeval- 
uated by a semi-empirical method 131 • In addition to the radiative correction to the weak matrix 
element, this method accounts for the related corrections to the energy levels and to the electric 
dipole (El) amplitude. Despite it is intended to incorporate the radiative and correlation effects, it 
is unclear how the results obtained by this method are related to those derived in the framework of 
the rigorous QED approach. The total binding QED correction obtained in Ref. amounts to 
-0.32(3)%. 

In the present paper we describe in detail the complete aZ-dependence evaluation of the one- 
loop QED corrections to the PNC transition amplitude in alkaline atoms and present the corre- 
sponding numerical results for the 6s-ls PNC amplitude in cesium [22] and for the 7s-8s PNC 
amplitude in francium, which is going to be a subject of the PNC experiment, as proposed in Ref. 
I24 I. Moreover, we perform a high-precision atomic structure calculation of the PNC transition 
amplitude in francium using a large-scale configuration-interaction Dirac-Fock (CI-DF) method 



26| ]. The obtained contributions are combined 



and compare the results with those from Refs. 
with other terms to improve the theoretical predictions for the PNC transition amplitudes in Cs 
and Fr. 

The relativistic units {h = c = 1) and the Heaviside charge unit (a = e^/ (47r), e < 0) are used 
throughout the paper. 



3 



II. QED CORRECTIONS 



A. Formulation 



A systematic derivation of the QED corrections in a fully relativistic approach requires the use 
of perturbation theory starting with a one-electron approximation in an effective local potential 



V{r) 



-iCt - V + Pm + \/(x))^/'„(x) = Enlpni^) 



(1) 



In neutral atoms, it is assumed that V{r) includes the interaction with the Coulomb field of the 
nucleus as well as partly the electron-electron interaction. The interaction of the electrons with the 
quantized electromagnetic field and the correlation effects are accounted for by the perturbation 
theory. In this way we obtain quantum electrodynamics in the Furry picture. 

To derive formal expressions for the transition amplitude we employ the method developed in 
Ref. I27'] and described in detail in Ref. Il28[l . While this method is valid for arbitrary many- 
electron atom and for arbitrary (single, degenerate, and quasidegenerate) initial and final states, its 
formulation is especially simple for a one-electron atom (or an atom with one electron over the 
closed shells) and for the case of single initial and final states. 

We consider the transition of the atom from the initial state a (which is 6s for Cs and 7s for Fr) 
to the final state b (which is 7s for Cs and 8s for Fr) accompanied by the absorption of a photon 
with momentum k, e nergy / c° = |k|, and polarization e*^ = (0, e). The transition amplitude is 



given by the formula ||27 



2811 



Zs'^^T^ i dE' I dEg,,^,a{E\E) 



X 



V2m 



dE gbbiE) 



-1/2 



1 



dE gaa{E) 



Ta 



-1/2 



(2) 



In the case under consideration (one electron over the closed shells), the Green functions 

{E', E), gaa{E), and gbb{E) are defined by 

gb-^,a{E', E)6{E' -E-oo) = J d^d^'^i{^')G^{E\ x'; 00; E, x)7Va(x) , (3) 

gaa{E)S{E' -E) = —j rfxf/xVl(x )G(i?', x'; E, x)7V,(x) , (4) 



gbbiE)6iE' - E) 



2ti 



dxdx.'^l{x.')G{E' , x; E, x)7Vfe(x) 



(5) 
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Here 

r\ 1 poo P 

GJE',x';uj;E,x) = — — — ^ / dx^dx'^ / d^y exp (iE'x'° - lEx^ - icuy^) 

xA^(y)(0|T^(x')j.(2/Mx)|0) (6) 

is the Fourier transform of the Green function describing the process, ^^(x) is the electron-positron 
field operator in the Heisenberg representation, ^ = ^/'^7°, 7° is the Dirac matrix, 

v/2A;0(27r)3 

is the wave function of the absorbed photon, and 

x'; x) = — — / dx°dx'°exp{iE'x'° -iEx°){0\Ti:{x')ip(x)\0) (8) 

is the Fourier transform of the Green function describing the atom. The contours and F^ sur- 
round the poles corresponding to the initial and final levels and keep outside all other singulari- 
ties of the Green functions. It is assumed that they are oriented anticlockwise. Z3 is the renor- 
malization constant for the photon wave function and the factors ^ dE gaa{E) and 

2^ dE gbb{E) serve as the normalization factors for the electron wave functions of the 
states a and b, respectively. The Green functions G and Gj are constructed by perturbation theory 
after the transition to the interaction representation and using Wick's theorem. The Feynman rules 



for G and G^ are given in Ref. 112 8h . 



To the lowest order, the PNC transition amplitude is described by diagrams presented in Fig[T] 
Denoting the contribution to gh.y a{E',E) from these diagrams by g^^}^^{E\E) and taking into 
account that g^} = 1/{E — Sa) and gf^^ = 1/{E — Eb) and, therefore, the normalization factors in 
formula are equal to 1 , we obtain 

= if / ^^CaiE'^E). (9) 

According to the Feynman rules [28] and definition we have 

{b\ea'^ A^\n) {n\Hw\o) 1 



2n E'-Sb^ E-Sn E-Sa 

n 

i I sr^ {h\Hu,\n) {n\ea'' Ay\a) 1 

2ti E'-Sb^ E'-Sn E-Sa' ^ ^ 

n 

Here 

Hw = -{GF/V8)QwPN{rh5 (11) 



is the nuclear spin-independent weak-interaction Hamiltonian 111, Gp is the Fermi constant, 75 
and a" = 7^7^ = (1, a) are the Dirac matrices, and pjq is the nuclear weak-charge density 
normalized to unity. Substituting expression (fTOI) into equation ® and taking into account that, 
for a non-Coulomb potential V{r), there is no states of different parity but of the same energy, we 
obtain 



-(0) 



E 



{h\ea.^ Ay\n) {n\Hw\o) ^ {h\Hu}\n){n\ea" Ay\a) 



(12) 



We note here that the case of degenerate levels, which takes place for the pure Coulomb potential, 
can be considered employing the related formulas from Ref. [28]. 

The one-loop SE corrections to the PNC transition amplitude are defined by diagrams presented 
in Fig. 121 To derive the formal expressions for these corrections, one has to expand formula Q to 
the next-to-leading order: 



-(1) 



dE' £ dEgi]lSE',E) 



2m 
1 

'2 



dE' f dEgllSE',E) 



1 

27n 



dE gil\E) + 



2m 



dEgl^^E) 



.(13) 



Let us consider the derivation of the contributions from the diagrams "a" and "c". According to 
the Feynman rules ll28ll . we have 

{b\i:{E')\ni){ni\ea''A,\n2){n2\Hw\a) 1 



£,Ae',e) 



271 E' - Eb 



{E' - en,){E - e. 



n2) 



E 



gi^''\E',E) = - 



— E 



{h\L{E')\n^)(rix\Ew\n2)in2\ed'A,\a) 



n\n2 



{E' - en,){E' ~ e. 



E-Sa 



(14) 
(15) 



Here the SE operator is defined as 



{c\nE)\d)^^j duj2 



{cn\I{uj)\nd) 

E — UJ — UEn 



(16) 



where I{uj) = e^a^a'^D^^(uj), D^^(uj) is the photon propagator defined as in Ref. ll28ll . and 
u = 1 — iO ensures the correct position of poles of the electron propagators with respect to the 
integration contour. Taking into account that, for a non-Coulomb potential, the energy in Eq. 
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([T4b is never equal to Ea and the energy Eq. (fTSt is never equal to Sb, we obtain 



f dE' f dEgi^;^liE',E) = -2m[j2 

''^b J'^a ni,n2 



{b\J]{eb)\ni){ni\ea''A^\n2){n2\Hw\a) 



(17) 



f dE' f dEgi;^fl{E',E) = -2m[Y, 



{b\T,{eh)\ni) {ni\Hw\n2) {n2\ea'' Ay\a) 



n2j 



n 

-E 



{b\i:'{eb)\b){b\Hw\n){n\ea''AM 
(6|S(£6) |6) {b\Hw\n) {n\ea''A^\a) 



(18) 



where T.'(E) = dT,(E)/dE. The contributions containing {b\Yl'{eb)\b) should be considered to- 
gether with the second term in equation dTSt . Taking into account that 



2m 



— i dEgil\E) = — i dE 



mmb) 



27ri Jr {E - SbY 



m{sb)\b) 



(19) 



we obtain 



dE' dEgflSE\E)\ [— ^ dEAgll\E) 



1 



(6|ea'^74,^|?T,) (n|ifiy|a) ^ (6|if^|n)(n|eQ;''y4,^|a) 



It, 

Adding this contribution to the terms dlTt and dTSt . we obtain 



(6|S'(£fe)|6). (20) 



r 



(l,a) 



?ii,n2 



(^b ^ni)(^a 



1 ^ (6|S'(efe)|6)(6|eaM^|n)(n|iyvK|a) 



(21) 



ni=ib 



(6| ^(ef,) (ni |'^2) ('"'2 lea^'A^ |a) 

(^6 ^ni)(^6 ^n2) 



1 >^ (&|S^(£fc)|&)(6|gvt/|ri)(ri|eaM^|a) 

'2^ Bb-Bn 

{b\i:{eb)\b){b\Hw\n){n\ea''A^\a) 



E 



(22) 
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Similar calculations yield 

^(1,6) 



= 



^(1,9) 



"^^'^ {b\Hw\ni){ni\ea''A,\n2){n2\J:{ea)\a) 



"2; 



1 (6|ifvF|'^)('^|eQ;'^Ai/|a)(a|S'(£a)|a) 
"^9 2^ 



2 £ft — £^ 

n 



_ {b\ea''A^\ni){ni\Hw\n2){n2\T,{£a)\a) 

ni,n2 

1 (6|eQ;'^yl,y|n)(n|i^M/|fl)(fl|S'(£a)|a) 

~'~2 ^ F — F 

n 

ea'^/l,^ I n) (n I i/vK I o) (a I E I a) 



-E 



(^a ^n)^ 



.(i,e) _ _ (6| eg'' I ni ) (ni | S (g^) | ^2 ) (na | if | a) 

(^a ^ni)(^a ^^12) 



.(1,/) _ _ v^" (6|ij"iy|ni)(ni|E(£fc)|n2)(n2|eQ;M^|a) 



ni,n2 



^ (6n2|/(c^)|nin) 
[£b - a; - M£ni][£a - a; - we^J ' 

^ r , ^ (6|//„-|//)(//i|a,'U.|n2) 



(nn2|/(a;)|nia) 



X 



[£6 - a; - li£ni][£a - 0^ - ti£n2] ' 



(6n2|/(ci;)|nin) 

[Sb-UJ - USniHSb -co - USn^l ' 

^ (nn2|/(a;)|nia) 
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r 



(i,fc) 



/oo 
-oo „_ . 



r 



{ni\ea''A^\n3){n3\Hw\n2) 

- - Mf^ngJ [£^a - - USriil 

iU) _ _ i f°° \^ {bn2\I{uj)\nia) 



/oo 
'OO „. „_ „ 



27r / „ fe;, — Co" — 1 

ni,n2,n3 ^ " 

(ni|ij';y|n3)(n3|eaM^|n2) 
X T TT T • ^-^ 

[£6 - - M5„3J - - M5„2j 

Taking into account the corresponding diagrams with the mass counterterm results in the replace- 
ment S(^) ^ Sr(E) = S(^) - 705m. 

Since the wave length of the absorbed photon is much larger than the atomic size, one can 
use the dipole approximation. It means the replacement exp (zk ■ x) 1 in the photon wave 



function and, therefore, ea^Aj^ |e| (a • e) / ^/2k^{27Ff in formulas dHl), (III])-<E21>- Within this 
approximation, the corresponding formulas in the length gauge are obtained by replacing ct with 
r in all vertices corresponding to the photon absorption and by multiplying the amplitude (fT2b 
with the factor i{Eb — Ea), where Ea = ea + {a\Y^-s^{ea)\(i) and E^ = eb+ {b\T.^{eb)\b), and the 
amplitudes (l2Tt - (l32b with the factor i{Eb — £a)- This prescription can be derived from equations 
(fT2ll . (|2ni - (l32b employing the commutation relation a = i[hj:,, r], where Hd = —ia. ■ V + /3m + 
V{r) is the Dirac Hamiltonian. Alternatively, one can get it using Eq. (205) of Ref. [28] and the 
equal-time commutation relations for the field operators in the Heisenberg representation. 

The theoretical and experimental results for the PNC amplitude in alkaline atoms are generally 
presented in terms of the -EpNc amplitude which is defined as the matrix element of the z com- 
ponent of the atomic electric-dipole moment between the initial (a) and final (b) s states with the 
angular momentum projections rua = m;, = 1/2. It is related to the r amplitude by the equation 

-EpNC = irlea^A^ ea;]/{Eb - Ea) = rlea^A^ -4] , (33) 

where the s states a and b have the angular momentum projections = rub = 1/2, Ea and Eb 
are their total energies, and = ez is the z component of the dipole moment operator (e < 0). To 
zeroth order, one easily finds 



^l{b\d;,\n){n\Hw\a) {b\Hw\n){n\dz\a) 
EPNC - —— + - 



(34) 
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The one-loop SE correction is given by the sum of the following terms: 



ni,n2 



(6| SR(£fe) |ni) (ni 14 [722) {n2\Hw\a) 



l^ (6|S-R(£,)|6)(6|4|n)(n|gvy|a) 



A??b _ ^'^"^ (6| if I ni) (ni I 1 712) (712 1 (gg) I g) 

PNC — / / \ / \ 

, lv- (6|gn l")("l''J")("|SBfa)l<') 



n 



(6|SR(£fe)|ni)(ni|ift^|n2)(n2|4|a) 



4E 



{h\T.'^{ek)\h){h\Hw\n){n\d,\a) 

{Sb — En) 

{h\L^{eb)\h){h\Hw\n){n\d,\a) 



XFd _ ^\^"N^|4|?T'l)(?T'l|^H/|"'2)("'2|SR(ea)|a) 
ni,n2 ^ '^-niA'^a '^"2.' 

1 ^ {h\dz\n){n\Hw\a){a\'^K{^a)\a) 
_ ^ (n|ifty|a)(a|SR(£n)|a) ^^^^ 



^^f _ (&|iiv^|ni)(ni|ER(£j,)|n2)(n2|4|Q) ^^^^ 



(5n2|/(a;)|nin) 

[£6 - a; - UeniliSa - co - UEn^l 
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27r 



{b\Hw\n){ni\d^\n2) 



X 



n,ni,n2 

{nn2\I{uj)\nia) 



(42) 



5EI 



PNC 



{ni\Hw\n2) {n\dz\a) 



X 



n,n\ ,n2 

{bn2\I{uj)\nin) 



[£b — ^ — UEn-i^] [^b — ^ — US. 



(43) 



"2 J 



5E{ 



PNC 



(6|4|^)(ni|-f^vy|n2) 



n,n\,n2 

{nn2\I{ijj)\nia) 



[Sa-UJ - USn^] [Sa-UJ - USr 



(44) 



^-^PNC 



(6r^2|/(u;)|^^la) 
[efo - u; - 



ni,n2,n3 

{ni\dz\n:i){n3\Hw\n2) 



PNC 



2^ 



{bn2\I{uj)\nia) 
[Sb-uj - n£.„J 



(45) 



ni,n2,n3 

(ni|i^iy|n3)(n3|4|n2) 



(46) 



According to Eq. ( B^ . the corresponding expressions in the velocity gauge are obtained by the 
replacement d^ —ieoLzj [Eb — Ea), where the energies Ea and Eb include the SE corrections. 
In addition to the replacement d^ — > —iea^/ (eb — £a) in Eqs. (I34b - (l46t . it yields the contribution 

(6|SR(£,)|6)-(a|SR(£j|a) 



r rpudd 
"-^PNC 



which results from the expansion 

1 1 



-E, 



£b — 



PNC : 



(47) 



Eb — Ea Eb — Sc 



{h\i:^{eb)\h) - {a\T.^{ea)\a) 



(48) 



{Sb - So) 

It can be shown that the sum of contributions (l43t - (l46t is the same in the length and the velocity 
gauge. Because of the gauge invariance of the total SE correction, the same is valid for the sum of 
the other terms, Eqs. (|35l)-(l42l) and (l47b . 
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Formulas (l35l) - (l47t contain ultraviolet and infrared divergences. To cancel the ultraviolet di- 
vergences, we expand contributions (l35t - (l40b into zero-, one-, and many-potential terms and con- 
tributions (HTt - iH^ into zero- and many -potential terms. The ultraviolet divergencies are present 
only in the zero- and one-potential terms. They are removed analytically by calculating these 
terms in the momentum space (for details, we refer to Refs. jlollBoljBlll ). For the standard zero- 
and one-potential terms we employ the equations given in Ref. f3(J\ whereas the corresponding 
expression for the zero-potential PNC term is presented in the Appendix. The many-potential 
terms are evaluated in configuration space employing the Wick rotation in the complex to plane. 
The infrared divergences, which occur in contributions (l35t - (l38t and (l45b - (l46b . are regularized by 
introducing a nonzero photon mass and cancelled analytically. 

The expressions for the VP corrections, which do not contain any insertions with the external 
photon line or the weak interaction attached to the electron loop, are obtained from Eqs. (I35t - (l40b 
by the replacement of the SE operator with the VP potential. The other VP corrections will not 
be considered here, since their contribution is negligible. To a high accuracy, the VP potential is 
determined by the Uehling term, which corresponds to first nonzero term in the expansion of the 
vacuum loop in powers of the Coulomb potential. The renormalized expression for the Uehling 
potential is 



oo 



Uvemir) = -az'^ j dr' Aixr' p{r') j dt{l + ^'^^ 

1 



2 



t2 

[exp (— 2m|r — r'\t) — exp (— 2m(r + r')t)] 



(49) 



Amrt 

where p(r) is the nuclear charge density, normalized to unity. To account for the screening effect 
on the Uehling potential, one should replace Zp{r) by Zp{r) — [Z — l)pcoTc{r'), where Pcorei^) 
is the charge density of the core electrons, normalized to unity. The higher-order one-loop VP 
potential, so-called Wichmann-KroU term, can be evaluated for the point-charge nucleus using 
approximate formulas derived in Ref. ll32ll . 



B. Local Dirac-Fock potential 

Since the energy intervals between the levels 6s, 6pi/2, 7s, and lpi/2 in Cs and the levels 7s, 
7pi/2, 8s, and 8pi/2 in Fr are very small, to get reliable results for the transition amplitudes under 
consideration, one needs to use a local potential V{r) that reproduces energies and wave functions 
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of these states on the Dirac-Fock (DF) accuracy level or better. We construct such a potential by 
inverting the radial Dirac equation with the radial wave function obtained by solving the DHF 



equation with the code of Ref. ll33ll . 



The radial DF equations have the form 113 3h 



\dr r J V r / 

i + ^\Ga+{Vc + ^)Fa-mFa =eaFa ^" 



r 



(50) 

ar r J \ r / r 

Here Ga/r = Qa and Fa/r = fa are the large and small radial components of the Dirac wave 
function of the a shell electron, Sa is the one-electron energy, = (— + 1/2) is the 
relativistic quantum number, Vq is the Coulomb potential induced by the nucleus, and Ya{r)/r 
is the screening potential. The functions and consist of two parts. The first part is the 
result of the action of the exchange-interaction operator on the radial wave functions Ga and Fa. 
The second part is the contribution from the non-diagonal Lagrangian multipliers, which provide 
the orthogonality of the radial wave functions corresponding to different values of the principal 
quantum number Ua but the same Ka. The functions X^ and X^ are calculated self-consistently 
from the DF equations employng the radial wave functions obtained at the previous iteration step. 
Let us consider the Dirac equation for the a shell electron with a local potential Va{r): 



(51) 



^ - -)Fa + Va{r)Ga + mGa = SaGa 

dr r J 

(I. + ^\Ga + Va{r)Fa - m Fa = BaFa. 

\dr r J 



In contrast to the nonrelativistic Schrodinger equation, generally speaking, it is impossible to 
choose such a local potential Va{r) which would exactly reproduce the one-electron energy £a and 
the radial components Ga and Fa for a given shell. This is due to the fact that the potential Va{r) 
enters both radial equations. However, one can derive an approximate potential by inverting the 
radial Dirac equation for the large component: 

vyr) . - ,„ ^ i- (1 - f„ . ^ 2^ . ^ A-f . (52) 

This leads to a local potential V^{r) which has some singularities, because the function Ga has 
nodes in the core region for > + 1. 

Let us consider another method of constructing the potential Va{r). Multiplying the first and 
second radial Dirac equations with Ga and Fa, respectivily, and summing them, we obtain 

-Ga (^-—]Fa+ Fa ( ^ + —] Ga + K(r) Pa + TrGI - mF^ = SaPa, (53) 
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dr r J \ dr r 



where pa = Gl + F^. Inverting this equation with respect to \4(r), we have 

T.(l)/ N Ga ( d Ka\ Fa f d Ka\ Fj 

Vy[r) = Sa ^ 3 }Fa -3- H }Ga + m m — 

Pa \dr r J pa \dr r J pa pa ^^^^ 

r Par 

Despite the potential Va^^ (r) has no singularities in the core region, it can oscillate and singularities 
can occur in the nonrelativistic limit. 

To smooth the potential Va^\r) in the core region, we use the following procedure. Instead of 
the density pa, we consider an average density p^ defined by 

PnaKa = ^nK^PuKa, ^ ^"'^^ ^ ^' ^^^^ 

n<na n<na 

where WnKa positive weights. Thus, the density Pua^a gets some admixture of the densities of 
the core shells corresponding to the same value of but different values of the principal quantum 
number n < Ua- Since the maximal values of the core shell densities are located nearby the nodes 
of the function Ga, the density can be made to be smooth and nodeless by a proper choice of 
the weights WnKa- Outside the core region the densities Pna^a ^iid p^^^a almost coincide with 
each other. This is due to a fast decrease of the core wave functions outside the core. Assuming 
the nonlocal part of the DF potential can be replaced by a local potential which is the same for all 
shells with the same one can derive 

VF\r) = Vc + ^^ + j\ Y (^«- + ■ (56) 

rriaKa n<na 

(2) 

The potential Va (r) derived for the shell a can also be used for all shells with the same and 
different values of Kq. This potential with weights WuKa oc (m — enaKa)/{^ ~ ^uko) was used in 
our calculations. 

In Table m we compare the energies of the cesium atom obtained with the local potential V{r), 
that was derived using mainly the DF wave function of the 6 s state, with the DF energies and with 
the experimental ones. The corresponding comparision for the francium atom, where the local 
potential was derived using mainly the DF wave function of the 7s state, is presented in Table HH 



C. Numerical evaluation of the QED corrections 



Numerical evaluation of expressions (l34l) - (l47t was performed by employing the dual-kinetic - 
balance finite basis set method ll42ll with basis functions construced from B-splines. The cal- 
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culation of the zeroth-order contribution (l34l) . with V{r) constructed as indicated above, yields 
EpNc = -1.002 for i33cs and Ep^^c = -10.19 for ^sspr, in units of ixlO-^^{-Qw) /N a.u. 
These values should be compared with the corresponding DHF values, —0.741 for ^'^^Cs and 

— 13.72 for 223pj.^ ^jj]^ jj^g values that include the correlation effects, —0.904 for ^'^'^Cs and 

— 15.72 for 223pj. ^ggg jj^g j^g^j section). The individual SE corrections are presented in Table 11111 
Since there is a significant cancellation between terms containing the infrared singularities, the 
terms corresponding to = a in Yl'^{ea) and n = b in T.'^(£a) are subtracted from contribu- 
tions (B5b - (BHl and added to contributions (H3l - (H3b . The total SE correction SE^^q, presented in 
Table Uni contains also the free term, — a/(27r)i?PNC5 mentioned above. Since this term is usu- 
ally included into the weak charge Qw, one has to consider the binding SE correction defined 
as 5Ep^Q = + 0!/{2n)Ep^c- According to our calculations, the binding SE correction 
amounts to —0.67% for cesium and —1.29% for francium. To estimate the uncertainty of these 
values due to correlation effects, we have also performed the calculations with V{r) constructed 
employng the DF wave function of the 7s state for cesium and the 8s state for francium. While 
this leads to a 2% decrease of the transition amplitude, the relative shift of the SE correction is, 
however, five times smaller. Since the correlation effects contribute to the transition amplitude 
on the 20% level, we assume a 4% uncertainty for the total SE correction. Therefore, our value 
for the binding SE correction is —0.67(3)% for cesium and —1.29(5)% for francium. In case of 

cesium, our value differs from the previous evaluations of the SE effect, which are —0.9(1)% ll20ll 
I — I 

and -0.85% [19]. 

We have also calculated the VP correction. The individual contributions for the Uehling part, 
calculated including the screening correction as described after equation (l49b . are presented in 
Table IIVI The total Uehling correction is almost independent of the screening effect and amounts 
to 0.410% for cesium and 1.037% for francium. These results agree well with the previous calcu- 
lations of this correction. The individual contributions for the Wichmann-KroU (WK) correction, 
obtained employing approximate formulas for the WK potential from Ref. I^^, are given in Table 
fVl The total WK correction is equal to —0.004% (cf. llA ) for cesium and —0.028% for francium. 
This leads to the 0.406% result for the total VP correction for cesium and to the 1.01% result 
for francium. Therefore, the total binding QED correction amounts to —0.27(3)% for cesium and 
-0.28(5)% for francium. 
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III. ELECTRON CORRELATION EFFECT ON THE PNC TRANSITION AMPLITUDE 



To calculate the correlation effects on the PNC amplitude we start with the relativistic Hamil- 
tonian in the no-pair approximation: 

H^^ = A+HA^, H = J2hDU) + Vl^' + Vi-\ (57) 

j 

where /id is the one-electron Dirac Hamiltonian, the index j = 1, ...N enumerates the electrons, 
and and are the Coulomb and the Breit electron-electron interaction operator, respec- 
tively. The frequency independent Breit interaction in the Coulomb gauge is given by 

i<j i<j 

(58) 

Here Vq^ is the so-called magnetic or Gaunt term and is the retardation term. The operator A+ 
is the projector on the positive-energy states, which is the product of the one-electron projectors 

A+ = A+(l)---A+(iV), (59) 

where 

-^+(0 = I M^)){un{i) I . (60) 

n 

Here Unii) are the positive-energy eigenstates of an effective one-particle Hamiltonian h^, 

KUn = BnUn, (61) 

which can be taken to be the Dirac Hamiltonian hr,, the Dirac Hamiltonian in an external field or 



the DF Hamiltonian in an external field 142 



44 



To calculate the -EpNc amplitude, we add the weak interaction to the full Hamiltonian: 

H{fi) = H + HwU) , (62) 

j 

where Hw is defined by Eq. (fTTT) . 

With the PNC interaction added to the one-electron DF Hamiltonian, one obtains the coupled 



equations, which are usually referred to as the PNC-HF equations 114611 . The linearization of these 
equations with respect to the parameter n would make them inhomogeneous. Since in our calcula- 
tions we do not perform such a linearization, the equations remain homogeneous. Li this case the 
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PNC amplitude can be calculated using the equation 

d 



-SpNC — 7^ 



(63) 

^l=o 



where D = J2i is the dipole moment operator and and are the many-electron wave 
functions of the initial and final states, respectively. They obey the equations 

= E.ifi) M/^/i) , Hi^i) vl>/(/i) = Eji^i) vl>/(;x) . (64) 

The many-electron wave functions \E'* and are represented by a large number of the configura- 
tion state functions (CSFs): 

a 

The CSFs ^•^^ are linear combinations of the Slater determinants, which are constructed from the 
one-electron wave functions Expansion (E5l) contains the CFSs of different parity, since 

the weak interaction is included in the Hamiltonian H{fi). 

The one-electron functions are obtained as eigenfunctions of the Dirac-Fock operator in 

the external field: 

It should be noted that the Dirac-Fock operator /idfI/w) depends on the parameter jj,, since the 
one-particle density matrix is constructed from occupied orbitals Un{n)- We can also consider the 
set of one-electron wave functions m° (/i) defined by equations 

hl{^i)ul{^i) = e'Mu'M , h'M = /idf(O) + fiHw, (67) 

where /idf(O) is the standard Dirac-Fock operator without the external field. 

The PNC amplitude can be calculated in the Hartree-Fock approximation by using only one 
CSF in expansion (l65t . Using equation (l63t and the wave functions uKfi), one obtains so-called 
Dirac-Fock value of the PNC amplitude. If the set of u„(yu) is used, the method, in principle, is 



equivalent to the PNC-HF method, which was used by different authors II47LI48LI4S 

In the large scale configuration interaction (CI) method the set of the CSFs for given quantum 
numbers JM is generated including all single, double, and the most significant part of triple 
excitations in the positive spectrum of the one-electron states In what follows, this method 

of evaluation of the PNC amplitude will be referred to as the PNC-CI method. 
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To obtain the set of the one-electron functions m„(;u) and we solve equations (l66b and 

(l67t using the finite basis approximation, 



Ur. 



a 

with the basis functions given in the central field approximation: 



jla 



<^a(r) = — I . (69) 

The representation (1^^ differs from the standard one by the factor iK This factor is introduced to 
make the one-electron matrix elements of the PNC Hamiltonian to be real: 



oo 

(a I tVtv \h) = - (-l)('<>-'»+i)/2 s^^_^^ s^^^^^ j dr [P, Q, - P,_ 



(70) 







where /a + + 1 is even. With this one-electron basis, the large scale PNC-CI matrix is also real 
and Hermitian. 

For the occupied atomic shells, the large Pa and small Qa components of the radial wave 
functions are obtained by solving the standard radial DF equations. For the vacant shells the 
Dirac-Fock-Sturm equations are used. For details of the Dirac-Fock-Sturm method we refer to 
Refs. HQ. The basis set containing the radial functions up to 17s, 16p, 12d, 7f, 5g, and 2h 
states was used in the calculations. 

In the calculations of the one-electron PNC matrix elements (ITOb we used the Fermi nuclear 
distribution 

where t = 2.3 fm. The parameters c and po were determined to reproduce the value of the nuclear 
mean-square radius Rn = (r^)^/^ and the normalization condition for pAf(r). 

In Table |VI| we present the results of our calculations of the PNC amplitude for Rb, Cs, and 
Fr. The results obtained by the DF method are given in the third column. Our DF value for the 
6s-7s PNC transition in Cs, —0.741, is in a good agreement with the values —0.742 HlOl] and 
—0.739 llli . which were obtained by the direct summation over the intermediate states. For the 
7s-8s PNC transition in Fr our DF value, —13.72, is also in a good agreement with the —13.56 
resuh obtained in Ref. |||2^. Our PNC-HF values, -0.138 for Rb and -0.926 for Cs, can be 
compared with the values —0.139 and —0.927, respectively, obtained by a similar method in Ref. 
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1 4911 . In the fifth column of the table, we present our PNC-CI values, which include the core- 
polarization correlation effects. The uncertainty of these values is estimated to be on the 1% level. 
For comparison, the most accurate results by other authors are listed in the sixth column of the 
table. In the second column we give the values of the nuclear-mean- square radius Rn, which were 
used in our calculations. They were obtained by the formula Rm = 0.836^^/3 + 0.570 [51]. In 



case of Fr, the corresponding results with Rn taken from Ref. 115211 are also presented. 

To calculate the contribution of the frequency independent Breit interaction (BI) to the PNC 
amplitude, we included the magnetic Vq^ and retardation terms in all stages of the calcula- 
tions. As the first step, the BI was included in the radial Dirac-Fock equations. We will refer this 
approach to as the Dirac-Fock-Breit (DFB) method. On the second stage, the BI was added to the 
Dirac-Fock-Sturm equations and to the Dirac-Fock Hamiltonian hu{fJ.) in the external field (1^ . 
This method of calculation of the PNC amplitude will be called as the PNC-HFB method. Finally, 
we added the BI to the many-electron Hamiltonian H(fi) in the external field and performed the 
large scale CI calculation. This approach will be called as the PNC-CIB method. To estimate the 
role of the retardation part of the Breit intaraction, we repeated all the calculations including only 
the magnetic (Gaunt) part of the BI and then took the difference with the PNC amplitude, 
which includes the total BI. 

In Table lVIll we present the magnetic Breit retardation Breit contibutions to 

the 6s-7s PNC amplitude in ^^^Cs and to the 7s-8s PNC amplitude in 223pj.^ obtained by different 
methods. The comparison of the total Breit correction to the PNC amplitude with the most accurate 
results by other authors are presented in Table IVIIII Finally, in case of francium, our PNC-CIB 
value amounts to -15.58(16) [Rn = 5.640fm] and -15.55(16) [Rn = 5.658fm] for ^sspr, and 
— 14.21(14) for ^^°Fr. They are in a fair agreement with the most accurate previous results ll26ll . 



-15.41(17) [Rn = 5.640fm] for ^^'"^Fr and -14.02(15) for ^i^Fr. 

IV. TOTAL PNC AMPLITUDES 

To get the total 6s-7s PNC transition amplitude in ^^^Cs, we combine the most accurate value 
that includes the correlation and Breit effects O, -0.902(5), with the -0.27(3)% binding QED 
correction, the —0.19(6)% neutron skin correction ll^, the —0.08% correction due to the renor- 
malization of Qiy from the atomic momentum transfer g ~ 30 MeV down to g = ligl . and the 
0.04% contribution from the electron-electron weak interaction 13. 0]- The analysis of accuracy 
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of the atomic structure PNC calculations is based on calculations of the hyperfine 

splitting, decay rates, and energy levels. As it was argued in Ref. IiqI . QED corrections to these 
uantities can be neglected on the 0.5% accuracy level. Using the experimental value for Ep^c/P 
and an average value from two most accurate measurements of the vector transition polariz- 
abilty, f3 = 26.99(5)a| flQHS, we obtain for the weak charge of ^^^Cs: 

Qw = -72.65(29)exp(36)th • (72) 

This value deviates from the SM prediction of —73.19(13) IstI by 1.1 a. 

In case of francium, combining our PNC-CIB values, —15.55(16) for 223pj. ^j^^ —14.21(14) 
for ™Fr, with the —0.28(5)% QED correction and the —0.08% correction due to the renormal- 
ization of Qw from the atomic momentum transfer g ~ 30 MeV down to g = IiqI . we obtain 
-15.49(16) for ^^^Fr and -14.16(14) for ^^opr. 

In summary, we have calculated the QED correction to the PNC transition amplitude in Cs and 
Fr. In addition, we have performed an independent high-precision calculation of the correlation 
and Breit interaction effects on the PNC amplitude in Fr. We have derived the weak charge of 
^^■^Cs, which deviates by 1.1 a from the SM prediction. Further improvement of atomic tests of 
the standard model can be achieved, from theoretical side, by more accurate calculations of the 
electron-correlation effects and, from experimental side, by more precise measurements of the 
PNC amplitude in cesium or other atomic systems, particularly, in francium 13, Si- 
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APPENDIX A: ZERO-POTENTIAL PNC VERTEX CONTRIBUTION 

The zero-potential PNC vertex contribution is defined as 

{b\Aw\a) ^J^sJ^s ^.(p')C(p',p)Vk(|p - Pl)^a(p) , (Al) 
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where p = {e, p) and p' = (e', p') are four vectors, 

r»,(p'.p) = -4™y p^7. _ „, 7°7' _ „, 7'p . (A2) 



oo 



q = |q|, r 



f Air 
Vw(q)=il / rfrpTvf?") exp (zq ■ r) = ?7 — / dr rpN(r) sin (qr) , (A3) 
i Q Jo 

|r|, and 77 = —{Gp/ V8)Qw- In equation (lAlll . it is implicit that e = and e' = Sb- 



One can easily express V^r{p' , p) in terms of the standard vertex function T"{p', p): 

rO.(p', p) = rO(p', p)75 - -[2£'m(Co + Cn) + 2emCu - m^j'Coh'' , (A4) 

71 

where the coefficients Co, Cn, and C12 are defined as in Ref. 1301 . After the isolation of the 
ultraviolet divergences in r°(p', p), the finite part is given by 

Ci?(p',P) = —[{A + 4m^Coho + i>'{B,e' + B2e)+i>{Cie' + C2e) + Di>'-foi> 

+H,e' + H^e - ^e'm{Co + Cn) - SemCiajT' , (A5) 



where all the coefficients are defined as in Ref. ll30l] . Integrating over the angles in Eq. dAll) . one 
can obtain 



(6|AH^,fi|a) = -^^i'^-^^5^^,,^J^^,^^ dp dp'py" J ^d^[Vw{q)Qi{p',P,OPkiO 
+Vw{q)Q2{p',P,0PLm, (A6) 
where Pi{^) is a Legendre polynomial, k = (— l)-'+'+^/^(j + 1/2), 

Qi = [A + 4m2Co + + 526) + £(Ci£' + Cae) + D^'e 

+ H2e - 8e'm{Co + Cn) - 8emCuMp')fa{p) 
+ [p'{B^e' + B^e) + Dp'e]Up')fa{p) 

+ [p(Ci£' + C2£) + Dpe']Up)Up) + Dp'pUp')Up) , (A7) 

Q2 = + 4m2Co + + B^e) + ^(Cie' + C2£) + De'e 

-Hie' - H^e + 8e'm{Co + Cn) + 8emCi2]Up')Up) 
+ [p\B,e' + 52£) + Dp'e]Up')Up) 

+ [p(Ci£' + C2£) + Dpe']Up')~Up) + Dp'pUp')Up) , (A8) 

^(p) and /(p) are the radial components of the Dirac wave function in the momentum representa- 
tion, defined as in Ref. iQ. 
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TABLE I: The binding energies of low-lying states in Cs, in a.u. Tlie experimental energies are taken from 
Ref. 



State Local potential DF Exp. 

6si/2 -0.13079 -0.12824 -0.14310 

6pi/2 -0.08696 -0.08582 -0.09217 

6p3/2 -0.08479 -0.08397 -0.08965 

7si/2 -0.05621 -0.05537 -0.05865 

7pi/2 -0.04251 -0.04209 -0.04393 

7p3/2 -0.04175 -0.04143 -0.04310 



TABLE II: The binding energies of low-lying states in Fr, in a.u. 

State Local potential DF Exp. [35. 36. 37. 38. 39. 40. 4j_1 

7si/2 -0.13640 -0.13271 -0.14967 

7pi/2 -0.08857 -0.08629 -0.09391 

7p3/2 -0.08199 -0.08071 -0.08623 

8si/2 -0.05740 -0.05626 -0.05976 

8pi/2 -0.04297 -0.04222 -0.04436 

8p3/2 -0.04071 -0.04023 -0.04188 
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TABLE ni: The SE corrections to the 6s-7s PNC amplitude in Cs and to the 7s-8s PNC amplitude in Fr, in 



%. The results are presented in both the length (L) and the velocity (V) gauge. 







Cs 
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0.15 




1.31 




1.11 


1.84 




1.35 




0.34 




0.40 


-0.36 




-0.23 




-0.38 




-0.32 


-0.64 




-0.51 




-1.29 




-1.53 


-1.21 




-1.46 


X IT'f 

OEp^Q 
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TABLE rV: The Uehling corrections to the 6s-7s PNC amplitude in Cs and to the 7s-8s PNC amplitude in 



Fr, in %. The results are presented in both the length (L) and the velocity (V) gauge. 







Cs 






Fr 


Contribution 


L-gauge 




V-gauge 


L-gauge 


V-gauge 


"^PNC 

<^-^PNC 
r iTiadd 

^-^PNC 
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n Q(Y> 

-0.107 
-0.188 
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1.037 


TABLE V: The Wichmann-KroU corrections to the 6s-7.s PNC amplitude in Cs and to the 7s-8s PNC 
ampUtude in Fr, in %. The results are presented in both the length (L) and the velocity (V) gauge. 
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TABLE VI: The E'pnc amplitude, in units of ixlO ^'^{—Qw)/^ a.u., calculated by different methods 
without the Breit correction. 



RNvm 



DF 



PNC-HF 



PNC-CI 



Others 



85Rb 5s ^ 6s 
i33Cs 6s ^ 7s 



223 



Fr 7s 



210 



Fr 7s 



8s 



8s 



4.246 
4.837 

5.640 

5.658 [52] 

5.539 
5.545 [52] 



-0.110 
-0.741 

-13.72 

-13.69 

-12.51 
-12.51 



-0.138 
-0.926 
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FIG. 1: Feynman diagrams for the lowest-order PNC transition amplitude. The wavy line terminated 
with a triangle indicates the absorbed photon. The dashed line terminated with a cross indicates the 
electron-nucleus weak interaction. 
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FIG. 2: Feynman diagrams for the SE corrections to the PNC transition amplitude. The wavy line ter- 
minated with a triangle indicates the absorbed photon. The dashed line terminated with a cross indicates 
the electron-nucleus weak interaction. 
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